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Abstract 


Scaling laws for possible outer planet magnetospheres 
are derived. These suggest that convection and its associated 
auroral effects will play a relatively smaller role than. at 
Earth, and that there is a possibility that they could have 
significant radiation belts of energetic trapped particles. 
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i) INTRODUCTION 

Of the outer planets, only Jupiter is known, from radio astronomical 
investigations, to have a magnetosphere, which requires that the pressure of 
the planetary magnetic field be sufficiently large to stand off the dynamic 
pressure of the solar wind flow. At present, observation does not rule out 
magnetospheric interactions with the solar wind with the other planets. Any 
well-conceived program for the exploration of the outer planets must therefore 
be prepared for the eventuality that one or more might have magnetospheres. 

This eventuality implies that the design of suitable complement of detectors 
for the exploration of the unknown magnetospheres must be considered. In 
the absence of hard experimental information, such considerations will rely, 
however unwisely, upon theoretical extrapolations upon what is known about 
Earth and Jupiter. In this spirit then, this paper presents a highly specula- 
tive discussion of hypothetical outer planet magnetospheres. We take what is-; 
reasonably well understood about the Earth's magnetosphere, what is guessed at 
about Jupiter's magnetosphere, and extrapolate to possible magnetospheres of 
Saturn, Uranus, and Neptune. The theoretician's 'point of view is adopted 
throughout. Since the attenuation of the solar wind with increasing heliocentric 
distance , implies that the magnetic moments of the outer planets need not be large 
for them to have magnetospheres, it does not seem unlikely a_ priori that they 
will do so. However, with the exception of Jupiter, their magnetic moments are 
completely unknown. Therefore, we will concentrate upon developing a set of 
relations which scale the outer planets' magnetospheres to their unknown magnetic 
moments and to the properties of the solar wind extrapolated theoretically to 
the appropriate heliocentric distance. In order to illustrate the implications 
of these scaling laws, we will then compute some properties of the outer planets' 
magnetospheres, based upon the assumption that their magnetic moments scale as 
their rotational angular momentum. Clearly this procedure looks only under the 
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larnppost where there is some light; yet the extrapolation of terrestrial physics 
is the only intellectual procedure available. Prudence dictates that we must 
expect it. to err. 

In Section 2, we scale the size of a magnetosphere to its planet's magnetic 
moment and heliocentric distance, assuming that the balance of forces at the 
boundary of the magnetosphere -- magnetopause -- is earthlike -- namely, a 
pressure balance between a vacuum dipole planetary field and the solar wind. 

We also estimate the strength of its internal convection flow assuming it is 
driven as at Earth, by magnetic reconnection at the nose of the magnetosphere. 

We discuss procedures by which the density of plasma of ionospheric origin trapped 
in the magnetosphere may be estimated. In Section 3, we discuss possible 
radiation belts of trapped energetic particles. Here the limitations of our 
method are most starkly delineated. It is a general truism about turbulent 
plasmas that they generate energetic particles in a variety of ways. Yet only 
one of the mechanisms suggested for the generation of the Earth ' s . radiation belts-- 
let alone the energetic particles in laboratory and astrophysical plasmas -- radial 
diffusion, can be scaled a priori to arbitrary magnetospheres. Therefore, we 
pursue the consequences of the only hypothesis we can make. In Section 4, we 
present with all humility a table of properties of possible outer planet magneto- ~ 
spheres, based upon the assumption that their dipole moments scale as their 
rotational momentum. At this point, our method of extrapolation of terrestrial 
physics leads us to a very illuminating contradiction, namely, that the effects 
of planetary rotation are likely to be much more pronounced at the outer planets 
than at Earth. This leads us to doubt, for example, that present calculations 
•of magnetospheric shape, and perhaps even scale size, are adequate for the outer 
planets, and to the speculation that the outer planets could have powerful 
radiation belts. These general conclusions may retain some validity even though 
our specific magnetic moment estimates may err greatly. 
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Limitations of space unfortunately force us to presume of the reader a 
reasonable working knowledge of basic magnetospheric physics. For general 
reference, however, we present in Figure 1 a schematic of the Earth's magneto- 
sphere in which various features to be discussed are put in geometrical per- 
spective -- the magnetopause, the boundary between the shocked solar wind 
and the magnetosphere, the bowshock standing upstream of the magnetopause, the 
plasmasphere where cold plasma of ionospheric origin corotates with the Earth, 
the Van Allen belts, and part of the geomagnetic tail. The view is of a slice 
through the noon-midnight magnetic meridian, and most of the geomagnetic tail, 
which is some thousand Earth radii long, is not shown. 
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2) SCALING OF EARTH-LIKE MAGNETOSPHERES 
2.1) Characteristics of the Solar Wind 
A planet P located at a distance r A.U. from the sun and within the helio- 
sphere boundary, has a magnetic moment Mp, radius Rp, and rotation period Tp. 
With this information, together with an appropriate scaling of solar wind 
parameters, we may outline a model of its magnetosphere, assuming only that 
it is Earthlike. We scale the solar wind number density N, flow speed u, and 
radial and azimuthal components of the solar wind magnetic field, and B Q 
respectively, according to standard theory (Parker, 1963), normalizing to values 
typically observed at r = 1 . 


u = constant - 4 x 10 cm/sec 


N = 7/r 2 cm ^ 


D 5y 5x10 Gauss 
B r = T = 


5y = 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 



R V /2 

B 0 3 


5y /l+r 2 - — (r >> 1) 

/ 4 r 

r 


(2.4) 


The thermal conduction of the solar wind beyond Earth is not well understood 
The simplest assumption, which may err, lets the electron and ion temperatures, 

T g and T^ separately scale adiabatically . 


T 

e 


7x 10 4o K 

4/3 

r 


2x 10 4o K 

i ~ 4/3 

r 


( 2 . 6 ) 


Scarf (1969) has discussed the expected characteristics of the. solar wind 
near Jupiter in more detail. In particular, he suggests that the temperature 
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anisotropy will reverse, so that near Jupiter the perpendicular temperature 
T a will exceed the temperature T n parallel to the magnetic field direction. 
As a consequence, different electromagnetic wave instabilities (Kennel and 
Petschek, 1966) than those encountered near Earth (Kennel and Scarf, 1968) 
would be expected to reduce the thermal anisotropies. 


2.2) Nose of the Magnetopause 

The nose radius D of the planetary magnetopause can be estimated assuming 
that the dipole field is essentially a vacuum field, whose moment is oriented 
more or less normal to the ecliptic plane. Then, according to Spreiter and Alksne, 
1969, the radial distance Dp to the magnetopause at the subsolar point is deter- 
mined by the balance of solar wind dynamic pressure and magnetic pressure, the 
dipole field having been doubled by magnetopause surface currents: 


M, 


2 1/6 


2ttpu 


(2.7) 


Normalizing to the Earth, and using (2.2) to scale the solar wind dynamic 
pressure, we find 





( 2 . 8 ) 


9 

where Dg = 10 Rg = 6.4 x 10 cm. The magnetohydrodynamic solutions for the shape 
of the magnetopause, which scale as the single parameter D, indicate that the 
distance between the local dawn and evening magnetopause, is 3D. The magneto- 
.spheric magnetic field at the nose of the magnetosphere is /Strpu^ - . 

The criterion D p = Rp defines the minimum planetary magnetic moment for 
which a magnetospheric interaction is expected. Mien Dp = Rp, the surface 
magnetic field pressure is just large enough to stand off the solar wind 
dynamic pressure. In units of the Earth’s magnetic moment. Mg, the minimum 
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planetary magnetic moment M p * is 


M * 
P 

M„ 



(2.9) 


Scarf (1969) has suggested that the gravitational interaction of the solar 
wind with the massive outer planets may significantly modify the flow configura- 
tion about the magnetosphere. We would expect significant modifications when the 
gravitational potential energy of an ion just beyond the bow shock exceeds 
its thermal energy, KT\, 


W G _ g P R P ’ , 

KT. (1.3)D p (KT./M.) 

where g p is the surface acceleration of gravity, and 1.3 D p is the expected 
distance to the bow shock (see Section 2.3). 


2.3) Characteristics of Planetary Bow Shocks 

Magnetohydrodynamic calculations (Spreiter and Alksne, 1969) indicate 
that a bow shock should stand a distance 0.3 D p upstream from the nose of the 
magnetosphere. Shocks are expected at all the outer planets since the Alfven 
Mach number remains constant and the sonic Mach number increases with r, based 
upon (2. 1-2. 7). However, the structure of the shocks encountered could differ 
from those at Earth. For example, a significant component of the Earth's 
bow shock is a large amplitude magnetic whistler mode wave train (see Fredricks, 
et al., 1970). In order to stand ahead of the shock in the solar wind, the 

whistler phase speed upstream must match the solar wind speed. Since the 
maximum whistler phase speed is ~ Ab /M g C^, where is the Alfven speed, 
whistler wave trains are possible when 
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C A < U 


i /mT/M 


r e 


J A. 


( 2 . 11 ) 


which is satisfied for r > 1 if it is satisfied at r = 1, since is indepen- 
dent of r from (2.1-2. 7). M./M is the ion to electron mass ratio. On the other 
hand, electron plasma oscillations, which do not play a role in the Earth's bow 
shock, could be important beyond r = 1 (Scarf, 1969). The minimum phase velocity 
of these waves is the order of the electron thermal speed a , where, from (2.8), 

o e 

1.4x10 

a e ~ " ~~ 2 / x ' cm / sec - Whenever, u/a > 1, electron plasma oscillations could 
r 0 

stand in the shock. Very little is known theoretically or experimentally about 
shocks which are supersonic to electrons. 


2.4) Reconnection on the Dayside Magnetopause 

Dissipative interactions leading to tangential stresses at the magnetopause 
are responsible for the geomagnetic tail (Axford, Petschek, and Siscoe, 1965; 
Dungey, 1961), the internal convection of plasma and magnetic field within the 
magnetosphere (Axford and Hines, 1961), and energetic particle bombardment of th 
auroral zone ionosphere by the convecting plasma (Kennel, 1969; Axford, 1969). 
Whether the dissipation is due to enhanced viscosity arising from plasma turbu- 
lence at the magnetopause (Axford, 1964), or to the resistive reconnection of 
solar wind field lines with magnetospheric field lines (Dungey, 1961; Levy, Pet- 
schek, and Siscoe, 1964), or both, has not been clearly established. However, 
it has been established that magnetospheric substorms (Akasofu, 1968), which are 
due in part to enhanced convective flow, do result from enhanced field-line 
reconnection, since they correlate with the solar wind field component anti- 
parallel to the Earth's dipole field (see Arnoldy, 1971, and the references 
therein). For this reason, we will evaluate only the consequences of reconnection 
and not turbulent viscosity. 

The electric field, imposed on the magnetosphere by reconnection, should be 
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proportional to where B. is the component of solar wind field anti-parallel 

to the magnetopause magnetic field. has considerable temporal variation, 
leading to temporally unsteady convection and substorms in the Earth's magneto- 
sphere. However, ought roughly to scale as the magnitude of the solar wind 

magnetic field. Assuming that the proportionality between the planetary con- 

U ^A 

vection electric field Ep and — — does not vary with heliocentric distance, we 
find 


•p— ~ 1/r (2.12) 

E 

where Eg - 1 kV/Rg is a typical terrestrial convection electric field. (2.12) 
scales identically as the estimate of Brice and Ioannidis (1970), who used a 
specific theoretical model of reconnection (Petschek, 1964) to scale Ep. 

We may estimate the solar wind energy input, W p , into the magnetosphere 
as follows. If b is the magnetosheath magnetic field downstream from the 
bow shock, then the flux of magnetic energy transported towards the magnetopause 

to be dissipated by reconnection into internal magnetospheric convection is 

cE 2 " 

roughly (^—) (b /8 tt) . The area of the dayside magnetopause is the order of 

2 

ttD , so that 


W r 


n 2 cE b 
ttD r — n — 
b 8r 


(2.13) 


When the bow shock is strong, b will scale as the solar wind field B, so that 
we may use (2.5), (2.8) and (2.12) to scale (2.13) 


W 


P-. 


W r 


m t 


2/3 


1 

r 


(2.14) 


17 

where Wg = 5 xTO ergs/sec (Axford, 1964). Since there exist no generally 
accepted theories or laboratory experiments which scale the reconnection rate to 
plasma parameters, the estimates (2.12) and (2.14) may err. However, it is 
dangerous to assume that no reconnection occurs at all. 
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2.5) Tail of the Magnetosphere 

Assuming Ep is approximately uniform, then the electric potential 4>p across 
the magnetosphere is approximately 3 EpDp, so that 

1/3 


*E 


m e ; 


•2/3 


(2.15) 


where <j>p - 30-100 kV is a reasonable value. Magnetic flux is transported into 
the magnetospheric tail at the rate F = c<{>p, where 


F = Cc<}> E ) 


( M p 1 


1/3 


r' 2/3 = 3 * 10 12 


M 


1/3 


-2/3 


Maxwells/sec (2.16). 


e 


where <j>g - 30 kV was chosen. Since magnetic flux cannot accumulate indefinitely, 
a second magnetic neutral line is expected in the magnetospheric tail, at 
which reconnection again occurs (Dungey, 1961; Axford, Petschek, and Siscoe, 
J.964.).. The reconnected flux wou.i.d then .he converted towards the nose of the 
magnetosphere to replenish that which has been stripped off, by reconnection at 
the nose, to feed the tail. In steady state, -the two flux transport rates must 
be equal. By analogy with the geomagnetic tail, a plasma sheet, containing 
energetic plasma heated by Joule dissipation during reconnection and other 
processes, would be expected planetward of the tail neutral line. Since there 
is presently no adequate understanding of the temperature and density observed 
in the Earth's plasma sheet, nothing concrete can be said about the density and 
temperature of any other plasma sheets, other than that the total plasma plus 
magnetic pressure must be constant across the plasma sheet. Energetic particles 
precipitating from the convecting plasma in the plasma sheet should produce 
aurorae in the high latitude planetary ionosphere (Kennel, 1969; Axford, 1969). 

It is not understood theoretically why the geomagnetic tail contains the 
flux it does; consequently, reliable estimates for the flux stored in other 



possible magnetosp'heric tails are not possible. However, let us suppose that 
the length of the tail Lp scales geometrically as the Earth's, which is some 
100 nose radii Dp long (Dungey, 1965), so that Lp = 100 Dp. An estimate for 
the steady state convection time T is then 


T 


c 


L p 100 D p 
u u 


and 


T C (P) 

T C (E) 




(2.17) 


where -T (E) - 1.5 x 10^ sec ~ 4 hours. 

During the time T reconnection would transfer a flux FT to the magnetospheric 

v C 

tail. In steady state, this is the flux stored in each lobe of the tail, 

Fp, so that • 


F„ 'v 

p 'V 


100 c4>D: 


u 


and using (2.8) and (2.16) 


r 

p 



% 




(2.19) 


We may estimate the tail magnetic field as follows. Beyond a distance 

Dp downstream from the planet, the magnetic field should be stretched out in 

„ . .... • 

a tail-like configuration, in two lobes, with field in the solar direction in 
one lobe and the antisolar direction in the other. The lobes are separated by 
the plasma sheet. Assuming the field is essentially a vacuum field, and 
therefore uniform across the tail cross-section, the magnetic field in the tail 
is then 

2F 

B - (2.20) 

ttR t ; 


where is the tail radius (2.20) corresponds to one of the basic assumptions 
in the flaring tail models of Tverskoy (1968) and Spreiter and Alksne (1969b), 
who assume the tail to be a cylinder bifurcated by a plasma sheet. Near the 
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In the terrestrial ionosphere, an auroral "oval" (Akasofu, 1964) of 
enhanced particle precipitation surrounds the polar cap. The area of the oval 
is comparable with that of the polar cap. A significant fraction of the. energy 
of the convective flow is dissipated as auroral precipitation and ionospheric 
heating. An upper limit for the energy input per unit area in the auroral oval 
u)p may be found by combining (2.14) with (2.20), whereupon 


■ 


[VI 

(I)-, i 

Is) 

E J 


1b e J 


,-2/3 


(2.25) 


2 

where - 1-10 ergs/cm -sec is a reasonable typical value. 

2.6) Corotation and Convection 

Figure (2A), from Brice and Ioannidis (1970), schematically illustrates the 
streamlines, in the magnetic equatorial plane, of the convective flow from the 
plasma sheet towards the nose of the Earth's magnetosphere. There are two 
distinct regions, of open and closed streamlines. The open streamlines are 
convective return of flux to the nose of the magnetosphere; near the Earth, 
where corotation dominates convection, the flow streamlines are closed. The plasma 
remains in this region long enough to approach thermal equilibrium with the 
ionosphere; the plasma density is consequently relatively high within the coro- 
tation region. Outside the corotation region, plasma escaping from the ionosphere 
is convected rapidly to the magnetopause where it is lost. Consequently the 
density is lower in the convection region (Brice, 1967; Nishida, 1967). The 
boundary separating the high density plasmasphere and the low density convection 
region is ordinarily quite sharp. Figure (2B) describes the calculated plasma- 
sphere at Jupiter, where corotation is much more powerful than at Earth. 

There are several means by which the relative importance of corotation and 
convection may be parametrized. For example, we. may compute the ratio of the 
convection time (2.17) to the rotation period Tp, 
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planet, the tail must join smoothly with the nose of the magnetosphere. Conse- 
quently, we take R^, - 1.5 D, the radius of the magnetopause on the dawn-dusk 
meridian. Thus, combining (2.18) and (2.20) we find 


B I b 


30 <f> 

uD 


2.25 x 


r ° 4 D 


( 2 . 21 ) 


b t (p) 

VE) 


1/r 


( 2 . 22 ) 


where B„(E) = 38y. The tail field should decrease monotonically with distance 
downstream approaching the value /SrrP^ at asymptotically large distances 
(Spreiter and Alksne, 1969b), where is the static pressure in the solar wind. 

We may now estimate the area A of the polar cap, the region of field lines directly 
connected to the solar wind, since the flux leaving (or entering) each polar cap 
must equal the flux in each lobe of the tail. Thus, 


A = 


100 c<f>D 
2uB„ S 


(2.23) 


where B^ is the equatoraal magnetic field at the sufface of the planet. 
2 

Assuming A - tiTq , where r^ is a characteristic dimension, then 


r 

x rv 


50c<f>D 

V 


0 ^ / n uB ^ » anc ^ ^ ie co -*- a ' t E tuc ^ e ^ of the boundary of the polar cap is roughly 


_ 1° / 50c<j>p 

~ / TT uB„ s 


(2.24) 


where r^/Rp has been taken small 



-13- 


t c (P)t e 

V E > T P 


1/3 


M r 


T 

1/3 E 


(2.26) 


T /T < 1 implies a dominant convection region, as for Earth, whereas T /T D > 1 
implies a dominant corotation region, as for Jupiter. Similarly, we may compare 
the magnitudes of the corotation and convection electric fields, E^ and E^ 
respectively. At the magnetic equator on a given txibe of force, we have 


2ir LR P B ^ 


„ R d B d 
2tt P P 


J CR 


T n . 2 
P cL 


(2.27) 


where L measures the distance in units of planetary radii. The minimum value 

M 

------ - ; - 2n P 

of occurs at the magnetopause, where L = D/R p and E Ep £ ^ : — j • The 

- 1 ' - ‘"' K J p cD 


ratio of the minimum corotation field to the convection field is, using 


(2. 12). and (2.8), 


“crW .. f Tp ^ 

- C rT £ ) 1 T P 'A » E J - 


(2.28) 


where A^ R (E) - 0.3-1. Again, when A Ep (P) > 1> corotation dominates. (2.28) and 
(2.26) are identical. 

Finally, the plasma energy density associated with rotational motion could 

distort the dipole field (Melrose, 1967; Brice and Ioannidis, 1970). This 

effect is measured by the corotation beta, 3^^, the corotation energy density 

divided by the magnetic energy density. For a given plasma number density 

p.. at the magnetopause, maximuizes in the dipole equatorial plane at the 
M LK 

magnetopause : 

Sir 1 


8 


CR 


B 


2 2 P M 


2ttD 

1 

'X, — 

’ P M ' 

2ttD ' 

T J 

n- 2 

l p s J 

uT p J 


(2.29) 


2 2 

where we have used B / 8 tt ^ pu . The ratio Q = 


2ttD 

uT„ 


scales as 
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2 
"O ■ 

2/3 

f^E) 

l m e j 

1 

• f P - 


(2.30) 


where 


10 


-4 


When 8^p approaches 1, we expect the magnetopause calculations 
referred to in Section 2.1 to fail because then centrifugal forces must be 

included in the stress balance. 


2.7) Plasma Density Profile 

Ioannidis and Brice (1971) have estimated the plasma density in the Jovian 

magnetosphere by a method which can be scaled to other planets. First, they 

noted that only photo-electrons have sufficient energy to escape over Jupiter's 

gravitational potential energy barrier. They then scaled the terrestrial flux 

2 

of photo-electrons deduced by Perkins and Yngvesson (1968) by a factor 1/r cos0, 
where 0 is the solar zenith angle at the foot of a given Line of force in the 

.ionusph.e±-.e . They ithen .'computed the -f Tux .and energy of escaping electrons, and 

* 

assumed that hydrogen ions would be pulled out of the ionosphere to ensure 
charge neutrality. From this, they could deduce a diffusive equilibrium density 
model, assuming no plasma loss from a given tube of force. 

This model predicts extremely large densities beyond L = 6. Therefore, it 
must be amended by the inclusion of loss processes, of which the most significant 
is the outward radial diffusion of cold plasma which is driven by interchange 

1 sets an 

upper limit for the plasma density, and in the absence of other loss mechanisms', 
determines the density. Figure (3) shows the plasma density profile computed 
by Brice and Ioannidis (1971) in this fashion. 

Scarf (private communication, 1972) has pointed out that this calculation 
can easily be extended to Saturn, provided Mp is known. Since Saturn's gravita- 
tional field and rotation period are comparable to Jupiter's, the flux of 


instabilities which set in when 8^p = 1. Thus the condition 8^,^ 


excaping photo-electrons will be 1/r - 1/4 as large at Saturn as at Jupiter. 
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Thus, near Saturn, the plasma density would be 1/4 that near Jupiter, and far 
from the planet would be determined by the condition B^ = 1. 
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3) RADIATION BELTS 

3.1) Radial Diffusion 

The origin of the energetic particles trapped in the Earth's magnetosphere 
is not completely understood quantitatively. Particles up to a few tens of 
KeV are injected during magnetospheric substorms, when rapid convection from 
the geomagnetic tail to the inner magnetosphere greatly compresses and heats 
the plasma (Axford, 1969). The maximum particle energy attainable by flow 
compression is given by the convection potential across the magnetosphere. 

Thus, from (2.15), convection should provide particles with energies not 
exceeding EgfM/Mgj^^r where E^ - 30-100 KeV. It is also likely that 

plasma turbulence can statistically accelerate particles to high energies 
within the magnetosphere (Kennel, 1969). Such mechanisms are poorly under- 
stood at present; anyhow, they cannot be scaled to other magnetospheres. It 
has •bear. ,s,ugg.c.s.t.c.d (.s.c.c Tv.sr.3J.v0y., .13.63., anu .Lire: rerereuCes Ln.erein) that the 
energetic component of the Earth's radiation belts is generated by injection 
of low energy particles in the outer magnetosphere followed by inward radial 
diffusion driven by variable electric and/or magnetic fields. If the field 
variations have sufficiently low frequency, the particles' first adiabatic 
invariant y = T 1 /B (where T 2 is the component of particle energy in motion 
perpendicular to the magnetic field) is conserved. Therefore, as particles 
diffuse from weak to strong magnetic field regions, their energy increases. 

If a typical magnetic moment can be estimated for particles injected at the 
magnetopause, then typical particle energies at any point in the dipole field 
can also be estimated from y-conservation. 

Since the ultimate source of the radiation belt particles is the solar 
wind, it is useful to compute the magnetic moment in the solar wind, based 
upon the flow energy density: 



(3.1) 
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y = ~ M^/B = 16r MeV/ Gauss 

If y is conserved for that small fraction of the impinging solar wind flux 
which not only traverses the shock and magnetosheath but penetrates the 
magnetopause boundary, we may use (3.1) to estimate the energy of radiation 
belt particles. On the other hand, should the particle’s magnetic moments be 
randomized by turbulence in the bow shock and magnetosheath, the magnetic 
moments of particles at the magnetopause could be somewhat smaller than (3.1). 
The maximum particle energy produced by radial diffusion will be of order 
yBp , where Bp is the planetary surface field. .If yBp > 0.5 MeV, electrons 
with sufficient energy to generate synchrotron radiation could be produced. 

The intensity of the radiation belts produced by radial diffusion is 
proportional to the fraction of the solar wind particle flux which diffuses 
across the magnetopause. How this occurs at Earth is not well understood., 
However, one thing is clear. The magnetopause of a rapidly rotating planet, 
will differ considerably from the Earth's. For example, if 3 CR = 1, the 
magnetopause could be subject to interchange motions, If there is rapid 
counterstreaming of corotating magnetospheric plasma and flowing magnetospheric 
plasma, then the growth of two-stream instabilities could increase the particle 
transfer rate. Clearly, the structure of the Earth's magnetopause cannot be 
extrapolated to the outer planets with confidence. 

The radial diffusion coefficient is determined by the power in time-varying 
electric and magnetic fields with periods comparable to the particle's azimuthal 
drift periods around the planet. .The particle drifts stem from three sources: 
convective electric field drifts from the combination of corotation and convec- 
tion, drifts due to gradients in the magnetic field strength, and drifts due 

to field line curvature (Hess, 1969). The time for a nonrelativistic particle 

* 

to drift once around the planet via the magnetic gradient drift is 
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2 2 
o L R d 
T , •> 2irc P 

t d ( i ° = — 


(3.2) 


where e = 5 * 10 ^esu, c = 3 * 10"*’*'* cm/sec, and LRp is the radial distance 
from the center of the planet to the particle. We assumed, for simplicity, 
that the particle has no velocity parallel to the magnetic field and so is ' 
confined to the magnetic equatorial plane. (3.2) may be suitably generalized 
to include parallel motions, and so the magnetic curvature drift, and to rela- 
tivistic particles (Lew, 1961). For particles with the same L, scales as 
2 

Rp /y , and so using (3.1), we find 


V L ; p ) ! 

T d (L;E) r 


(3.3) 


where Tp(L;E) 


- 0.15 L hours. Since is proportional to L for a given y. 


the radial: diffusion .coefficient when magnetic -drifts .predominate depends upon 

different frequency components of the fluctuating electric and magnetic fields 

at different L; similarly particles with different y resonate with different 

frequency components at a given L. • 

A number of mechanisms to drive radial diffusion have been suggested. For 

example, low frequency variations in the convection electric field, due to 

a variable solar wind, are thought to be important for the Earth (Falthammar, 1965; 

Birmingham, 1969; Cornwall, 1971; Mozer, 1971) and have been considered for 

Jupiter by White (1971). For a given y and L, the diffusion coefficient D is 

2 2 , 

of order CE (u^)/B (L) , 'where Efw^) is the electric field amplitude at the 

drift frequency and B(L) is the equatorial magnetic field strength. At low 

frequencies, the electric field amplitude should be reasonably uniform spatially: 

if furthermore the frequency spectrum is reasonably smooth, then D - in a 

-3 

dipole field, since B(L) 'v L . Perturbations of the magnetospheric magnetic 
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field stemming from irregular magnetopause motions, again driven by solar 
wind variations have been considered for Earth by Nakada and Mead (1965) 
and for Jupiter by Chang and Davis (1962) and Hess and Mead (1971). This 
mechanism has a basic L"^ dependence. Consequently, both electric and 
magnetic diffusion tend to be weak on the inner L-shells, where the highest • 
energy particles are involved. As mentioned earlier, interchange instabilities 
driven by corotation have been suggested by Ioannidis and Brice (1971). 
Interchange instabilities are one member of the class of low frequency drift 
instabilities which could drive radial diffusion (Kennel, 1969; Cornwall, 1970). 
No specific diffusion rate or L-dependence can be estimated for these mechanisms, 
since they depend upon knowledge of the nonlinear saturation levels of the 
instabilities. 

When corotation dominates the magnetic drifts, so that T CR /Tp << 1, 
the energetic particles circle the .planet in approximately one rotation period. 

In this case, electric and magnetic field amplitudes at the corotation 
frequency determine the diffusion coefficient for particles over a wide range 
of both y and L. It seems much more likely that the -time-varying fields will 
not stem from irregular solar wind variations, but will be relatively more 
coherently driven by corotation itself. It stands to reason that radial diffu- 
sion could be quite efficient in corotation-dominated magnetospheres . 

One such radial diffusion mechanism has been proposed by Brice (1971) and 
Brice and McDonough (1972) for Jupiter. Solar illumination creates a periodic 
heating of the planetary atmosphere which leads to tidal wind systems. The 
winds then couple to ions in the dynamo region of the ionosphere to drive Hall 

currents; polarization of the Hall currents then leads to electric fields which 

map along magnetic field lines out into space. The net electric potential 

WB p SR p 

associated with the atmospheric dynamo is of order — , where W, a typical 

wind velocity is of order one-tenth the sound speed. For Jupiter, this potential 
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is roughly 10 MV, so fluctuating fields greatly- in excess of that expected from 
solar-wind irregular convection may be possible. Furthermore, the diffusion 
coefficient may have a much weaker L-dependence than those associated with 
solar wind variability. For both these reasons radial diffusion near the 
planet, could turn out to be surprisingly efficient in co-rotation dominated 
magnetospheres . 


3.2) Loss Mechanisms 

High frequency fluctuations near the particles' cyclotron frequencies, which 
Violate the magnetic moment invariant, act as a loss mechanism by slowly 
diffusing the particles in pitch angle until the magnetic moment is sufficiently 
reduced that they are not reflected by the dipole field gradients and so are 
lost to the atmosphere. An upper limit to the stably trapped particle fluxes 


U ^-t, JT— 

<A 41 J.gJl J-.A Vk.[UOJiV.) 


instabilities. One such limit, involving electromagnetic ion cyclotron and 


whistler instabilities has been calculated for the Earth's radiation belts 


(Kennel and Petschek, 1966; Cornwall, 1966). When this fluxes of electrons and 
protons trapped in the radiation belts are sufficiently intense, the anisotropy 
in their velocity distributions, which is maintained by loss of small pitch 
angle particles to the atmosphere, permits whistler and ion cyclotron waves to 
grow unstably. The resonant particles driving the instabilities are diffused in 
pitch angle by the waves to the atmospheric loss cones, whereupon they are lost 
to the atmosphere. Thus, since the instability reduces the fluxes of trapped 
particles, an upper limit for stably trapped fluxes is given by the flux Jp* 
required for marginal stability (Kennel, 1971) 
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J* is independent of the particle mass (and is consequently identical for 

electrons and ions) and of the background plasma density N. However, only 

electrons and ions which can resonantly interact with the waves can be diffused; 

this condition implies that only particles with energies greater than the 

2 

magnetic energy per ion pair, B /8irN, at the magnetic equator, will be scattered 
by whistler oz - ion cyclotron waves. 

Pitch angle diffusion can only reduce particle fluxes to the stably trapped 
limit when the precipitation lifetime is less than characteristic radial diffu- 
sion time. The minimum precipitation lifetime (Kennel, 1S69) , a lower limit to 

4 

the precipitation lifetime, scales as L and so is large on distant L-shells. 

Thorne and Coroniti (1971) have arrived at an upper limit model for the intensity 

of the Jovian radiation belt. They assumed electric field diffusion, of 

.crrp-ngth -to permit .part-id es to .diffuse past Tc, and that ir» j cct’Son 

at the magnetopause was sufficient to create particle flux above the stability 

threshold for whistler and ion cyclotron waves. Beyond L = 6-8 radial diffusion 

is faster than precipitation, and the particle fluxes to the stably trapp d 

limit in the range L = 6-8. Thus, instabilities serve as a valve limiting the 

injection of particles to the inner L-shells in this model. Near the planet, 

2 

B /8trN exceeds the expected particle energy yB, using the Ioannidis and Brice 
(1971) plasma density model, so that whistler and ion cyclotron waves may be 
stable near the planet. However, there remains the possibility that electro- 
static instabilities of the loss cone type (Rosenbluth, 1965) could act as 
a turbulent loss mechanism. Such instabilities, with frequencies appropriate 
to scatter electrons, have recently been discovered in the Earth's radiation 
belts (Kennel et al, 1969), but as yet our knowledge of them is insufficient 
to permit extrapolation. 
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3.3) Role of Satellites in Magnetospheric Physics 

Planetary satellites whose orbits lie within the magnetosphere, offer the 
possibility of studying a completely new type of flow interaction. For example, 
Io, which is located at L = 6, well within the Jovian magnetopause, presents 
an obstacle which is supersonic to the corotation flow. It should have an 
interesting interaction and wake, whose properties are barely guessed at 
presently. Hess and Mead (1971) have argued that Io can absorb radiation belt 
particles as they radially diffuse across Io's orbit; for the highly L-dependent 
and consequently weak radial diffusion coefficients applicable in the drift- 
dominated regime, it is likely that particle absorption by Io would drastiaally 
reduce the particle fluxes. However, this leaves the problem of accounting for 
the observed synchrotron radiation created by relativistic electrons within the 
orbit of Io by some other mechanism. If radial diffusion is to provide for 
synchrotron radiation, then Io cannot be a barrier. Recently, Hess and Mead 
have estimated the radial diffusion coefficient from the synchrotron emission 
profiles (Berge, 1966); upon extrapolating this to the orbit of I o, they 
find a significant fraction of the radial diffusion flux can get past Io (Hess, 
private communication). 

Io produces decametric radio emissions which have been discussed by Goldreich 
and Lynden-Bell (1969). They argue that the v x B electric potential, from 
corotation of the plasma, across Io's diameter, is the order of 0.5 MV. This 
large potential then drives magnetic field aligned currents in the tubes of 
force intersecting Io, which close in the Jovian ionosphere. These field-aligned 
currents then produce instabilities in the Jovian ionosphere, creating waves 
.with frequencies up to the electron cyclotron frequency at the foot of the field 
line. Similar waves, also apparently associated with field-aligned currents, 
have been observed in the Earth's ionosphere as auroral hiss. Since the dissi- 
pation of field-aligned currents may heat the ionospheric plasma and also produce 
energetic beams of "runaway" electrons. There is the interesting possibility 
that Io could be a source of radiation belt electrons. 
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4) HYPOTHETICAL MAGNETOSPHERES 

Table I lists the orbital radius r (in astronomical units), the planetary 

radius R p , the ratio R p /Rg, the planetary rotation period T p , the ratio T p /T p , 

the planetary magnetic moment M p in units of the Earth's magnetic moment M p , and 

s 

the surface field B p for the planets Earth, Jupiter, Saturn, Uranus, and 

s 

Neptune. All the parameters but M p (and B p ) are well known, and have been taken 
from Newburn and Gulkis (1971) . The Jovian magnetic moment Mj has been estimated 
from radio astronomical evidence (Warwick, 1970). It is not known whether Saturn, 
Uranus, and Neptune have magnetic moments. However, current understanding of 
the dynamo theory of planetary magnetism indicates it would be dangerous to 
presume they have no magnetic moment, since they are rapidly spinning objects, 
with a reasonable possibility of having a conducting liquid core. Furthermore, 
the minimum magnetic moment M* (Eq. 2.10) for which a magnetospheric interac- 

_ i , 1 * . 

•'bicm -will ; occur is ■ -quite small, M '* - TO Mg, -My* - 8 x 10' " M £ , and M^* * 5 x 10 
M p . For the purposes of illustration, we have scaled the planetary magnetic 
moments according to the "Magnetic Bode's Law" (Moroz, 1967) whereby the magnetic 
moment is proportional to the total planetary angular momentum, a rule which works 
fairly well for Earth and Jupiter. These estimates of M p greatly exceed the 
minimum moment required for a magnetospheric interaction. We have not performed 
any scalings for Pluto, since it is sufficiently small that it may not have a 
magnetic moment. 

Table II lists basic parameters defining the magnetospheric configuration: 
the nose radius D p normalized to the Earth's nose radius D p and also to the plane- 
tary radius R p , the length of the geomagnetic tail L p in astronomical units, and 
B^, the magnetic field strength at the nose of the magnetosphere. These hypo- 
thetical outer planet magnetospheres are much larger than the Earth's, both in 
absolute units and in units of planetary radii. The estimated length of Jupiter's 
magnetic tail is significant on the solar system scale. Should the nose radius 
estimates be correct, then the satellites JV, Io, Europa, Ganymede and Callisto lie 
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within Jupiter's magnetosphere ; Janus, Mimas, Enceladus, Tethys, Dione, Rhea, 
Titan, and Hyperion within Saturn's; and Triton within Neptune's; all of 
Uranus' satellites lie within its magnetosphere. A rich variety of satellite 
interactions with planetary magnetopsheres may therefore exist. 

However large or small the magnetic moments, and consequently the magneto- 
spheres of the outer planets may be, the magnetic field in their outer regions 
will be considerably weaker than at Earth, due to the attenuation of the solar 
wind dynamic pressure with increasing heliocentric distnace. The estimate of 
B^j allows us to infer that Neptune's surface field need exceed only 1 or 2y 
for it to have a magnetospheric interaction with the solar wind. 

Table II also lists parameters defining internal convection: the electric 

potential tf>p across the magnetosphere; the net energy input Wp from the solar 
wind into the magnetosphere; and the energy flux Wp into the high latitude 
ionosphere from auroral dissipation of the convective flow. Jupiter and 
Saturn should have considerably larger convection potentials, and absorb consi- 
derably more energy from the solar wind, than Earth, whereas Uranus and Neptune 
are comparable to Earth as far as convection is 'concerned. The auroral 
particle energy fluxes into the high latitude Jovian ionosphere could consi- 
derably exceed those at Earth. 

Table III lists parameters necessary for the comparison of corotation and 


convection: the convection time T , the ratio of convection to corotation time 

c 

T /T „ D , which also yields the relative ratio of corotation to convection electric 

C LK 


fields at the magnetopause, A^p, and Qp which characterizes 


6 cr at 


the magnetopause. 


Both A^p and Qp favor corotation at the outer planets relative to Earth. On this 


basis, then,. we expect relatively large regions of corotation flow, and rela- 


tively small regions of convection in these magnetospheres. Furthermore, Q p 
is more than an order of magnitude larger for all the outer planets than for 


Earth, which suggests the strong possibility of corotation- induced distortions 
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in the magnetic field and/or interchange instabilities at these planets. This 
indicates that the simple calculation of the nose radius based upon an undis- 
torted dipole field is incorrect and can at best be regarded as an order of 
magnitude estimate. Finally, the magnetopauses of the outer planets could be 
irregular and noisy, thereby permitting injection of more particles into the 
radiation belts than at Earth. 

Table IV lists several parameters of interest in radiation belt physics: 

the characteristic magnetic moment p in the solar wind, the maximum particle 

energy attainable by radial diffusion pBp , in MeV; the maximum particle 

energy attainable by convection, e<j>„; the drift .time T„ (L) in hours, the ratio 

T d /T cr of drift to corotation time, and J*L^ where J* is the stably trapped flux 

limit defined by (3.4). If radial diffusion could bring particles to the surface 

of the planet without loss, it would produce several hundred MeV particles at the 

outer planets. Only protons would achieve such high energies, since electrons 

would lose energy to synchrotron radiation. 'Hie rings of Saturn should sweep out 

any radiation belt particles, so the maximum particle energy expected at Saturn 

is probably an order of magnitude smaller than the 600 MeV listed. The ratios 

T d /T cr listed in Table IV indicate that beyond L = 2, Jupiter's radial diffusion 

should be corotation dominated; Saturn's, beyond L 3; Uranus, beyond L = 10, 

and Neptune's, beyond L = 14. The Earth's radial diffusion, by comparison, is 

never corotation dominated. When the corotation domination region extends close 

to the planet, as for Jupiter and Saturn, our previous arguments lead us to 

suspect that there may be efficient generation of high energy particles. 

Finally, the stably trapped flux limits are fortuitously similar for all the 

planets. Whether or not the stably trapped flux limit applies depends upon 

2 

whether the particle energies exceed B /8 ttN -- which implies a knowledge of the 
plasma density N -- and whether the minimum precipitation lifetime is less than 
the radial diffusion time. 
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There will undoubtedly be many surprises in the magnetosphere of Uranus; the 
rotation axis is inclined roughly 98° to the normal to its orbital plane, so that 
twice per orbit, its rotation axis points nearly towards the sun. This con- 
figuration will occur in 1988. This suggests the possibility of a new and unusual 
magnetospheric configuration, if it turns out that its magnetic moment is aligned 
more or less along its rotation axis, as is the case for Earth and Jupiter. 

W.P. Olson has calculated the shape of the nose of the Uranian magnetosphere based 
upon this assumption; his results are presented in Figure (4) . In this case, 
the "polar cap" points directly towards the sun, and there exists the possibility 
of direct penetration of solar wind to the planetary surface. Whether or not this 
implies an especially intense radiation belt is unclear. G.L. Siscoe (1971) 
has discussed convection and the topology of a possible Uranian magnetic tail: 
his results are presented in Figure (5) . Magnetopause and tail reconnection 

,hnt,’n .take n.l ace -on lines of force connecting to the magnetic pole. It seems "**' 

♦ 

likely that the atmospheric tidal dynamo, postulated by Brice and McDonough to 
drive radial diffusion at Jupiter, will be most unusual. at Uranus. 
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5) DISCUSSION 

We have scaled the magnetospheres of the outer planets according to the 
theoretical variation of solar wind parameters and to Moroz's magnetic 
Bode's law for the magnetic moments. In the absence of better information, 
we have assumed where necessary that Earthlike physics prevails at the outer 
planets. Yet this procedure suggests that the magnetospheres of the outer 
planets could be very different from the Earth's. Several broad conclusions 
emerge. First, convection and its associated auroral precipitation should 
play a relatively smaller role at the outer planets than Earth. Corotation 
dominates. This in turn suggests that solar wind particles may penetrate a 
disturbed magnetopause and radially diffuse into the dipole more efficiently 
than at the Earth. The general increase in solar wind magnetic moment with 
increasing heliocentric distance indicates that the radiation belt particles 
could be considerably more energetic than at Earth. Thus, the outer planets' 
magnetospheres could be radiation-belt dominated. The outer planets offer the 
possibility of studying satellite-magnetospheric Interactions ; the interaction 
of Io with the Jovian magnetosphere is known to be significant. Uranus may 
have an unusual magnetic topology, since its magnetic dipole axis is currently 
directed more or less in the solar wind directions. These hypothetical outer 
planet magnetospheres seem to 1 be more like each other than the Earth. In parti 
cular, a careful study of the magnetospheric processes at Jupiter may lend 
greater insight for the rest of the outer planets than further extrapolation of 
Earth-like physics. The above parametric studies suggest that the investiga- 
tion of the magnetospheres of the outer planets can be carried out using 
instrumentation already developed. 
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FIGURE CAPTIONS 


Figure 1. EARTH'S MAGNETOSPHERE 

Shown here is a slice through the noon-midnight meridian of the Earth's 
magnetosphere, With the relative geometrical locations of various features to 
be discussed subsequently in the text. 

Figure 2. CONVECTION AND COROTATION AT EARTH AND JUPITER 

(A) is a schematic of the streamlines of the Earth's convective flow in 
the magnetic ' equatorial plane, taken from Brice and Ioannidis (1970). Local 
magnetic times are indicated, with the solar direction, local noon, at the 
top of the figure. The region of closed streamlines, the corotation or 
plasmashere, contains relatively dense cold plasma of ionospheric origin. 

(B) is a similar schematic for Jupiter. 


Figure 3. PLASMA DENSITY IN THE MAGNETOSPHERE OF JUPITER 

This figure is taken from Ioannidis and Brice (1971). The dashed line 
approaching infinity near L = 10 is the result of loss-less diffusive 
equilibrium calculations; the dotted line indicates the density limit set by 
recombination, and the dashed lines labelled Bgj = 1, 10, 20g indicates the 
density limit set by interchange instability for various values of Jupiter's 
surface magnetic field B^j. Centrifugal effects confine these densities 
largely to the Jovian magnetic equatorial plane. 


Figure 4. NOSE OF THE MAGNETOSPHERE OF URANUS 

Shown here are results of calculations by W.P. Olson of the nose of 
Uranus' magnetosphere. The solar wind impinges ■ upon the planet at 0°. Magnetosheath 
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plasma could directly penetrate the magnetosphere along the 0° line. 

Figure 5. CONVECTION AND THE MAGNETIC TAIL IN THE MAGNETOSPHERE OF URANUS 
Reprodiiced here is a schematic of the convective motions postualed by 
Siscoe (1971) . The numbers label a magnetic tube of force successive instants 
in its interaction with Uranus. Point 2 corresponds to field annilation at 
the nose of the magnetosphere; N.S. denotes neutral sheet. Corotation around 
the dipole axis has been neglected; it would be expected to give the field 
lines a helical twist. 
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